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ABSTRACT 


Under the assumption of radiative equilibrium in the core, of 
convective equilibrium in the outer parts of the envelope, and of a pro* 
nounced inhomogeneity in composition, models of the interior of the 
Sun are derived. As sources of energy, the proton-proton reaction and 
the carbon-nitrogen cycle are considered. The derivation of sun models 
is carried through in three steps. First, a homogeneous model is de- 
termined for the initial state of the Sun, when the hydrogen reaction is 
just beginning. Second, on the basis of this initial model, the transmu- 
tation rates and the present composition are computed for every point 
in the Sun. Finally, using this inhomogeneous composition, a model 
for the Sun in its present state is constructed. 
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Symbol 

T 
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M(r) 

L(r) 

K 

€ 

a 

y 

a 

c 

H 

k 

G 
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Definition 

temperature 

total pressure 

density of matter 

mass of a sphere of radius r 

radiant power (energy crossing a sphere of radius r per 
second ) 

mass absorption coefficient (area per unit mass) 

power generated per unit mass 

energy liberated by thermonuclear reaction 

ratio of specific heats 

radiation density constant 

velocity of light 

mass of unit atomic weight 

mean molecular weight of matter 

Boltzmann constant 

gravitational constant 

polytropic index 



I. INTRODUCTION 


Mass M, radius R, and luminosity L of the Sun are considered 
as known. The distributions of density p, pressure P, and temper- 
ature T in the interior are to be derived. 

The assumption is made that the chemical composition of the 
solar material may be specified by the abundance of hydrogen and helium, 
X and Y, respectively, and by the abundance of the heavy elements, Z, 
which is related to X and Y as shown in the equation; 

Z = 1-X-Y. 

The sources of energy generation in the Sun are thermonuclear 
reactions. 

In the temperature range between 10^ and 10^ °K, the main 
source of stellar energy is provided by the proton-proton (pp) ch.ain, the 
carbon-nitrogen (CN) cycle, and the helium reactions (3o'), 

Below 10^ the pp chain is dominant; above 2X10^ °K, the 

CN cycle dominates; in the intermediate range, both reactions compete 
with each other. The 3a reactions come into operation near 10** °K, 
and therefore do. not apply to the Sun, The pp chain occurs in the fol- 
lowing stages: 


(a) 

,H' + 


+ V 

(b) 



^He^ + y 

I (c-1) 

He^ + He’ 
2 2 

+ 2 

II(c-2) 

,He’ 

1 .He- 

= ^Be'^ + 7 

(d-2) 

..Be’ 

4 


i- ^ V 

(e-2) 



: 2^He" 


Ill (c-3) 


(d-3) = gB® + 7 

(e-^3) 5B® ^Be®''*‘ + + 1/ 

(f-3) 4Be®"‘'-> 2 lie* . 

The total amount of energy Q available per pp chain is: 


Q„_= 26. 207 MeV 
PP 


1 MeV = 1 . 6021 X 10"® erg. 




i { epi 


The energy released per gram of hydrogen^estroye^ per 
gram of helium formed equals: 

q"’* = . 6 . 3 X 10^® erg /gram. 

PP 4H 

The CN cycle occurs in the following steps: 


<a) 

,C‘2+ ,H> 

7 

(b) 


-f jS + 1 / 

(c) 


= N‘®-f 7 
7 

(d) 


- 0 ^® + 7 
8 

(e) 

q 15^ ]v^l5 

gW ^ IN 

+ 0 ^ + V 

(f) 


= sC*" + ,He 


The result of the cycle is the conversion of four protons into a 
helium nucleus; the carbon nucleus with which the cycle starts is recovered 
unchanged. The total amount of energy obtained from the cycle is: 

= 25. 026 MeV. 


The energy released per gram of hydrogen destroyed or per 

gram of helium formed equals: 

, Q 

= 6.0 X 10^® erg/gram. 

The physical picture of a spherically symmetric star in equili- 
brium is expressed mathematically by four simultaneous, nonlinear, 
ordinary differential equations of the first order, which represent the 
radial gradients of P, M(r), T, and L(r). Since there are more than 
four unknowns, additional relations between the unknowns are deter- 
mined from the physical properties of the. material. These relations 
are called the "constitutive equations, " which involve usually the chemi 
cai composition. The four fundamental equations of stellar structure 
are: 


(1) the equation of hydrostatic equilibrium: 

dP ^ GM(r)P 

dr * 

which balances the gravitational force by the pressure gradient; 

(2) the equation of conservation of mass: 

=4.r*p, 

dr ^ ■ 

which expresses the mass of a spherical shell of radius r and thick- 
ness dr; 

(3) the luminosity equation: 

dh(r) 


dr 


= 47t rzpe , 


which expresses the radiant energy of a spherical shell of radius r 
and thickness dr; 

the temperature gradient equations: 
dT i L(r) 


dr 

dr 


4ac 


T3 47t r2 


T-i 


P dr 


for radiative transport, 


for convective transport, 


which express the law of conservation of energy. 


The constitutive equations comprise the equation of state, the 
equation of opacity, and the equation of energy generation. 

Since everywhere in the interior of the Sun the radiative pressure 
is 80 small in comparison with the gas pressure that it can be neglected, 
the equation of state for a perfect gas applies to a sufficient degree of 
approximation: 


P 


h P 
H 


(j = l/[2X + -i-Y + 2Z], 


It is not necessary to consider the degeneracy of matter, 

II. DERIVATION OF SUN MODELS 
In the beginning, the pp chain will be considered as the source 
of energy generation, and the assumption will be made that the chemical 
composition is uniform throughout the Sun, This is true only when it 
ust begins to generate energy by thermonuclear reactions. As the Sun 
ages, thermonuclear reactions produce chemical inhomogeneity between 
the core and the envelope. At first, a homogeneous composite model 
will be considered. An inhomogeneous model will be discussed later. 

A . Homogeneous Composite Model 

Solar models incorporating the proton-proton reaction and com- 
puted under the assumption of radiative equilibrium give ^ computed 
luminosity much higher than the observed luminosity. The reason for 
this discrepancy is that a hydrogen convection zone extends ir|to the 
interior. 

The central temperature is smaller in a star with an outer con- 
vection zone than in a star with no such convection zone. In a convective 
region the temperature gradient is lower than the radiative temperature 
gradient. An oujber convection zone will reduce the central temperature 
and luminosity. 


At the surface, the atmosphere of the S\in is in radiative cquiii< 
brium. At some depth, the radiative temperature gradient becomes 
unstable against convection. 

The radiative core equations are: 
dP _ GM(r)P 


KP L(r) 

T» 4irr* 


dM(r) 


= 4 Trr^ 


dUr) 


= 4rr^€(r) 


P = p T 

pH ^ 


K = 6 , 52 XiO*^ (z+ T"*’* 

€ = 2. 8X 10’” X* p T^*, 

The convective envelope equations are: 

dP _ _ G M(r)p 
dr r* 

dM(r) , 
dr 

= 4»r»p€(r) (; 

dr ' 

P = — — p T 

PH 

y 

P = KT ^ = KT**. 

In order to reduce the equations to dimensionless form, the 
Seiiwar/.schild variables p, t, q, f, and x will be introduced: 


G M* 

^ 4 7T 

pH GM 
k R 


M(r) = q M 

L(r) = fL 
r = xR, 


S 


This transformation reduces the equations (2) to the form; 

^ = pq 


dx 

ia. 


tx' 


px^ 


dx t 

p = Et^’* 

E = 47t K 


/2Lf* 

Ik/ 




(4) 


and the equations (1) to the form: 


. pq 

dx tx^ 

dt 

dx 

r 

t*‘2* 



df 

= D p* X* t 


dx t 

dx 


(5) 

3X6.52X10^< (z 

C “ ^ 

59. 3/ 


L 


4ac 

1 

Q 


„ _ 2. 8X10-MX* 


M*-’ 4., 


4 7T 

\ k / 

LRW ** . . 



Since the structure of the outermost parts of the Sun is uncertain, 
the model of the envelope is varied. A series of models is computed 
corresponding to a series of values of E, Finally, the value of E is 
fixed by the condition that the overall model of the Sun fits the observa- 
tional data, 

»* • 

To describe the structure of the outer convective zone, the 
<.;f|uaiions (4) arc integrated with the boundary conditions: 
at,x 1: q = 1, p = 0. 

Eight homogeneous models were constructed, covering the range 
E - 1. 68 to E = 44. 821# Since the models refer to the homogeneous initial 

i 

state while the observed luminosity and radius refer to the present state 
of the Sun, no reliable results will be obtained in this way. These models 


are applied to the observed mass, luminosity, the radius of the Sun by 

selecting a value for the hydrogen content X, and then determining the 

values of ihe helium content Y, and the convection parameter E, with 

the help of the mass -luminosity relation and the energy-output relation. 

The main result obtained is the fact that a value for Z of 0. 024, 

in agreement with the spectroscopic data, is reached when E has the 

small value of 0, 86, that is, for a convective envelope extending only 

over 10. 9 percent of the radius. The completely radiative models yield 

much smaller values for the abundance of the heavier elements than the 

one w’hich is observed spectroscopically. 

The Sun is in radiative equilibrium at its center. It does not 

have a convective core. Therefore the equations of radiative equilibrium 

may be integrated from the center outward to the point at which they 

become convectively unstable according to the criterion: 

- / d log P \ ^ / d log P 

\ d log T \ d 


= 1 dlog T /rad” ^ d log T )adiab “ ^’adiab’ 


From this point outwa-rd the structure is convective. 

The boundary conditions at the center are: 
at X = 0: q =: 0, f = 0. 

The parameters C and D can be removed from the equations (5) 
by the following transformations: ; 


P 

q 

f 

t 


Po P 

% 'l’ 

f r 

0 

‘o 


X = 


X X 
0 
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In terms of these variables, the equations (5) assume the form: 




*>v 


-Sr-^ 

"•X 


t'’' 


dr 

dx" 

dr 

dx'" 


1,1 Sf\ 








X 


'•2 '>'2 A . 

p X t 


i\iZ.5 


% . 


to ^ 

p« 

^0 0 
^0 ^0 

c 


^ 1 


p 1.75 f 
^0 0 


f- 9'25v 
^0 ^0 


= 1 


D ~r = 1. 


( 6 ) 


If the value of one of the five unknowns , f ^ , t^ , x^ is 

selected arbitrarily, the remaining four can be determined by the con- 
ditions in the equations (6) in terms of the chosen parameter and of C 
and D. The central value of t is chosen: 


t — t 

‘'0 ‘'C* 


In addition to the variables x, p, t, q* and f, the following 
functions are computed: 


U 


V 
W = 
n + 1 = 
n + 1 = 


d log M (r) _ 4 TTi^ p 

d log r M(r) 

- d log P _ p GM(r) 
d log r 


P r 

^ ^96 4„^sp±llL. 

dlogr Ur) 

d log P _ Km ac GM(r) 
d log T " 3 P(r)/C (r)L(r) equilibrium, 

7 


for radiative 


7+ 1 


= 2, 5 for convective equilibrium , 


(7) 


The actual numerical integration of the equations (6) is performec 
by employing logarithmic variables: 


log p" = 

X 

log r = 


log q'’’ = 


log x'*' = y . 

(8) 

log r =: 

T 
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In terms of these variables, the equations (6) are reduced to the 

form: 

log(-dX/dy) - {jj - T ~y 
log (dip/ dy) = X-^-T+3y 

(9) 

log {-dr/ dy) = 1.75X - 9. 25 t + ^- y 


log {d(j)/ dy) = 2X + 2.5T - 0h3y, 


The combination of the equations (7), (8), and (9) yields: 

U = d(Jj/dy 

V n? -dX/dy 

(10) 

W = d0/dy 

n+1 = [dX/dy ] / [dr/dy] , 

The equations (9) are integrated under the boundary conditions 
at the center: 


at x’" = 0: q"' = 0, f- = 0, r' = 1, p"' = p* 
for a number of values of p^, which is treated as a parameter » The 
values of U, V, W, and n+1 are tabulated together with those of q'*% 

T'% and p'*' against the values of x'‘S starting from x'*' = 0. The out- 
ward integrations are stopped at x'*' = xj' , where n+1 attains the value 
2. 5. In this way the U - V curves are obtained along which n + 1 attains 


the value 2. 5 for the radiative solutions. This procedure yields the 


value of the inner solution at the interface p*‘^, q-’i,. t-i,. f'"' 


1C* 


ic 


lie 


1C* 


1C* 


^ic point where the convective envelope has to be fitted. 


Then the equations (4) are integrated under the boundary conditions: 


at X = 1 : t = 0, q = 1 , 

for a number of values of the parameter E, The values of p, q, t, U, 
and V are tabulated against x. In general, no convective envelope U - V 
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curve passes through the end point ^ic radiative core U - V 

curve. Then an interpolation is performed in the U - V plane between 

the convective outer solutions to find the value of E for the solution 

which passes through the point ^ic* procedure determines 

the envelope solution uniquely. For this value of E, the values of the 

outer solution at the interface p. 

X s 


IS 


X 


IS 


With the values known for x 


, t. , q. , 

and X. 

are found from the 

' is' ^is' 

IS 


= u. , V. 

= V-. 

Therefore: 

1C' is 

ic 


• t- = t 
ic» Hs ^0 

^i’c' ^is 

>.1^ 

~ % ^ic • 


- ' ^ic* Pic' ^is' ^ic' ^is' 
q?^, the values x^ , p^ , t^ , and q^ are determined. The equation for 

df/dx in equations (5) is then integrated for the envelope, using the solu- 
tion for the value of E which fits it at the end point of the inner 

solution for the assumed value of p'^. Since at 
X = 1, f = 1, and at X = x. f = f. : 

AO AS 

1 - f^g = D Jp2t2-Sx^dx 
^is 

f. = 1 - D f p^t^-Sx^dx, 

IT 

This equation determines f^^ at the interface for a value of D, 
which is determined by the following procedure. The continuity of f 

at the interface yields; 

_ 

f. = f. = f fr 

IS AC 0 AC 


^0 ^ic ~ ^ ^ / P^t2*5x2 dx . 

^^is 


By solving this cqAiation and the condition equation for D in equa- 
tions (6), with the knowledge of values of , p^ , , f'" , the values of 

xc« 
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f and D are obtained; 
0 


1/ [f. 


1C p2^2.5jj^3 
^0 0 0 


J dx ] 


X. 

IS 


° = ^0 /tPj 


Thci equation for C in equations (5) yields the value of C. Terms 
p and p"' yield p , and t = t gives the central value of t. 

The method of fitting the solutions indicates that a definite initial 
value of p^ determines C, D, and E uniquely. On the other hand, the 
fitting of the radiative core at the point beneath a convective 

envelope with a prescribed value of E determines C, D, and p^ uniquely. 
The term E is, therefore, the parameter which determines the model 
uniquely. 

B. Inhomogeneous Composite Model 

The initial homogeneous model of the Sun consists of a radiative 
core and a convective envelope whose depth depends on the parameter E. 
Energy is generated in the radiative core where the material remains 
unmixed. The rate of energy generation i's proportional to the density 
and to some power of temperature. It decreases with increasing dis- 
tance from the center, so that the hydrogen abundance varies at different 
rates at different points on a radius inside the core. 

Since at the center the transmutation rate from hydrogen to helium 
is about ten times higher than the average rate, a pronounced inhomo- 
geneity in composition exists at present in the Sun. Therefore, an inho- 
mogeneous model has to be constructed. 

The core equations are the same as the equations (1). For a 

chemically homogeneous envelope, the equations are the same as equations 

$ 

(2) except that /I has to be replaced by For a given value of E, the 
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envelope equations can be integrated under the boundary conditions: 
atx=l: q=l, t = 0, 
to obtain a unique solution. 

The energy and opacity laws assume the form: 


€ = € U 

oe f*e k 




K tf <|0*75/ ii \°'^%75_0.7S T”**- 
«e f'e I k M P ^ 


-4.25 


2. 8X10"” 


K,.e= 6 , 52 X 10 ’“" (z + -i^) (1 +Xg) 


0.75 


( 11 ) 


i = < 


1 = i 


m 

1 + X 
1 + X, 


i = 


JL, 




The equations (3) are replaced by the transformation: 


M(r) = qM 

^ P 47T R4 

L(r) = fL 


r = X R 

P =J2- JL 


t l^e 


R 

M 

4;r R’ 


(12) 


When equations (12) are substituted in equations (1)', the following 
replacement is obtained for the equation (5): 

dp/dx - t (pq/tx*) <3t/dx = -Cij<’-”(p*”r/t''-“ x^) 

dq/dx = i (px^/t) df/dx = Diip^t^’^x^ 

1.25 


D 


^^oc 

1 

4ac 

(4;r)^-” 

^oe 

/ GH r-' 

47T 

\ k / 


.7.5 


L.R 




(13) 


M 


6.5 


LR 


7.5 


M 


4.5 
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In terms of the starred variables, whose values have to be 
taken frona the core solution of the initial model, the equations (13) 
become: 

dp*/dx* = -i (p*q*/t*x**) 
dq*/dx’' = t /t* 

drVdx" = - i 

dfVdx’'‘ = i ip*' t*'-‘ X*'). 


( 14 ) 


with the boundary conditions: 

at x'*' = 0: q'*'= 0, f *' = 0, t'** = 1, p'*' = p'‘\ 

The condition equations are the same as in the equations (6). 

In the solution of equations (14), p^ is a trial parameter, as in 
the case of the equations (6) for a homogeneous model. Because of the 
presence of i, j, and I in equations (14), has to be considered as an 
additional parameter. The value of q^ may be taken from the homoge- 
neous model, and then it may be improved by integrations. 

The envelope solution forms a family of curves corresponding 
to the parameter E, The inner and outer solutions may be fitted in the 
U - V plane, as in the integration of the homogeneous model. 

In order that the effect of the carbon-nitrogen cycle upon the 
structure of the interior may be taken into account, the energy genera- 
tion is .represented by different laws, the differential equation governing 
the energy generation is modified, the non-dimensional parameter D is 
redefined, and a quantity 6 - 1 representing the ratio of energy generated 
by the carbon cycle to that generated by the proton-proton reaction 
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is introduced: 


€pp = 9.5X10-»OX^ p 
e CN = 3 XlO->« T 20 

df/dx = ij2Dp^t^x^6 


6 - 1 


D 



6, 3 X Z 

3 X9. 5 X10"^°X^ 


/ HGM .. 16 

l~k^/ 


9.55 X10“^ 

47T 



hlSJ 




The integration from the center is performed by the same pro- 
cedure which was explained previously. A 3 in the case of inhomogeneous 
composition, in addition to (n + 1)^,, becomes a second \mknown para- 
meter, The equation for 6-1 also contains the unknown quantity tg , 

By the use of the condition equations in (6), t^ may be eliminated and 
6-1 expressfed in terms of , M, R, X^, and Z. (Refer to Table I. ) 
The radiative solutions obtained by the integrations are fitted to the 

convective solutions which were previously obtained. 

* < 

TABLE I. SOLAR PARAMETERS 


Mq = 1.989 XI o3°kg 

A0 = 4. 5 XIO^ years ' 

Rq = 6. 9598 X 10® m 

X = 0. 74 

L^ 3.90 X 10^^ watt 

0 

Y = 0.24 

p 0 = 1. 409 g/cm^ 

Z = 0. 02 

T^ 0 = 5800 "K 

/X = 0. 599 

M, . _= 4.72 
bol © 

J/ 

i 

y spectral classification dG 2 
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Ill, RESULTS AND CONCLUSIONS 
The physical properties of the models obtained are summarized 
in Table II, from which the following conclusions can be drawn. There 
is a large abvindance of elements heavier than hydrogen and helium 

t ■« 

(Table III). The value lies between 1 and 10, The radiative core 

contains about 99 percent of the solar mass and extends over 84 percent 
of its radius. The convective envelope does not penetrate more than 
16 to 20 percent of the radiu/ 4 . The central density of the present Sun 
is 136 grams/cm^; its central temperature is 14X10^°K. 


TABLE U. ABUNDANCE OF ELEMENTS IN THE SUN 


Atomic 


Log of Relative 

Atomic 


Log of Relative 

No. 

Element 

Abundance 

No.’ 

Element 

Abundance 

L 

H 

12. 00 

38 

Sr 

2. 80 

2 

He 

11. 16 

39 

Y 

3.01 

3 

Li 

0. 86 

40 

Zr 

2. 16 

4 

Be 

2. 35 

41 

Nb 

2. 00 

6 

C 

8. 56 

42 

Mo 

1.96 

7 

N 

7.98 

44 

Ru 

1. 36 

8 

O 

9.00 

45 

Rh 

0. 70 

9 

F 

5.4 

46 

Pd 

1.00 

10 

Ne 

8.44 

47 

Ag 

-0. 50 

11 

Na 

6. 30 

48 

Cd 

1.76 

12 

Mg 

7.28 

49 

In 

0. 80 

13 

A1 

6.21 

50 

Sn 

1.2 

14 

Si 

7.60 

51 

Sb 

1.82 

15 

P 

5. 44 

56 

Ba 

2. 26 

16 

S 

7. 17 

57 . 

La 

1.8 

17 

Cl 

5.4 

58 

Ce 

2.4 

18 

Ar 

6.62 

59 

Pr 

0.6 

19 

K 

4. 96 

60 

Nd 

2.0 

20 

Ca 

6. 38 

62 

Sm 

1.5 

21 

Sc 

3. 20 

63 

Eu 

1.4 

22 

Ti 

4.96 

64 

Gd 

1.1 

23 

V 

4. 03 

66 

Dy 

1.6 

24 

Cr 

6. 00 

68 

Er 

0. 1 

25; 

Mn 

5. 30 

’ 69 

Tm 

0.5 

26 

Fe 

6. 76 

70 

Yb 

1.42 

27 

Co 

4.74 

71 

• Lu 

1.0 

28 

Ni 

5. 80 

72 

Hf 

0.4 

29 

Cu 

4. 99 

73 

Ta 

0.0 

30 

Zn 

4.53 

74 

W 

0. 2 

31 

Ga 

2. 16 

76 

Os 

0. 5 

32 

Ge 

3. 31 

77 

Ir 

-0. 2 

37 

Rb 

2. 28 

78 

Pt 

1.6 




82 

Pb 

2.68 


15 


TABLE m. PHYSICAL PROPERTIES OF THE SUN 


X 

Y 

Z 

E 

X. 

1 

^i 

(°K X 10^) 

Pi 

(°K XlO^) 

Pq 

0. 80 

0. 19 

0. 01 

2. 05 

0.849 

0.998 

0.94 

0.018 

14. 1 

1 36 

0. 74 

0. 24 

0. 02 

2. 14 

0. 847 

0.998 

1. 00 

0.019 

14. 6 

134 

0. 70 

27 

0. 03- 

2. 19 

0.846 

0.998 

1. 04 

0. 020 

15. 1 

132 

0. 60 

0. 32 

0. 08 

2. 33 

0. 842 

0.998 

1. 17 

0. 023 

16. 4 

128 


It has been found that most of the energy comes from the region 
located within 0 . 25 of the center of the Sun, The energy is transfer- 
red to the surface of the Sun by radiation and convection. Radiation due 

t 

to X- and 7 -rays is the dominant process out to about 0, 7 Rq» with con- 
vection becoming the principal process from 0. 7 to 1. 0 R^ . 
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